Phase-space structures in quantum-plasma wave turbulence 
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Abstract 

The quasilinear theory of the Wigner-Poisson system in one spatial dimension is examined. Con- 
servation laws and properties of the stationary solutions are determined. Quantum effects are 
shown to manifest themselves in transient periodic oscillations of the averaged Wigner function 
in velocity space. The quantum quasilinear theory is checked against numerical simulations of 
the bump-on-tail and the two-stream instabilities. The predicted wavelength of the oscillations in 
velocity space agrees well with the numerical results. 
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I. INTRODUCTION 



Quantum plasmas have attracted a renewed attention in recent years. The inclusion of 
quantum terms in the plasma fluid equations - such as quantum diffraction effects, modified 



equations of state 



and spin degrees of freedom [3] - leads to a variety of new physical 
phenomena. Recent advances include linear and nonlinear quantum ion-acoustic waves in 
a dense magnetized electron-positron-ion plasma |4], the formation of vortices in quantum 
plasmas [5|, the quantum Weibel and filamentation instabilities jg]- 10], the structure of weak 
shocks in quantum plasmas 11], the nonlinear theory of a quantum diode in a dense quantum 
magnetoplasma 12], quantum ion-acoustic waves in single- walled carbon nanotubes l^ . 
the many-electron dynamics in nanometric devices such as quantum wells Q, Q, the 
parametric study of nonlinear electrostatic waves in two-dimensional quantum dusty plasmas 



16j, stimulated scattering instabilities of electromagnetic waves in an ultracold quantum 



plasma [17], and the propagation of waves and instabilities in quantum plasmas with spin 
and magnetization effects 18l. Il9|. 

However, to date, only few works have investigated the important question of quan- 
tum plasma turbulence. The most notable exception is the paper by Shaikh and Shukla 
2o| . where simulations of the two- and three-dimensional coupled Schrodinger and Pois- 
son equations - with parameters representative of the next-generation laser-solid interaction 
experiments, as well as of dense astrophysical objects - were carried out. In that work, 
new aspects of the dual cascade in two-dimensional electron plasma wave turbulence at 
nanometric scales were identified. Nevertheless, the quantum-plasma wave turbulence re- 
mains a largely unexplored field of research. A reasonable strategy to attack these problems 
would consist in extending well-known techniques issued from the theory of classical plasma 
turbulence in order to include quantum effects. 

In this context, the simplest approach is given by the weak turbulence kinetic equations 
first derived in Refs. 2l|]-24|, the so-called quasilinear theory. In quasilinear theory, the 
non-oscillating part of the distribution function is flattened in the resonant region of velocity 
space. It is interesting to carry over the basic techniques of the classical quasilinear theory to 
the kinetic models of quantum plasmas. The resulting quantum quasilinear equations would 
be a useful tool for the study of quantum plasma (weak) turbulence, and for quantum plas- 
mas in general. In this regard, it is natural to initially restrict the analysis to the quasilinear 
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relaxation of a one-dimensional quantum plasma in the electrostatic approximation. 

Some earlier works explored the similarities between the classical plasma and the 
quantum-mechanical treatment of a radiation field 25J-29]. In those papers, one of the 
aims was to obtain information on the classical plasma through a quantum-mechanical 
language. For instance, the relaxation of an instability can be viewed as the spontaneous 
emission of "quanta" of a radiation field described by some quasilinear-type equations. How- 
ever, the application of the quasilinear method to a truly quantum plasma governed by the 
Wigner-Poisson system (i.e. the quantum analogue of the Vlasov-Poisson system) seems to 
be restricted to the work of Vedenov [301]. Surprisingly, there has been no systematic analysis 
of the consequences of the quasilinear theory in the Wigner-Poisson case. The present work 
is a first attempt in this direction. 

This manuscript is organized in the following fashion. In Sect. 2, the quantum quasilinear 
equations are derived from the Wigner-Poisson system. In comparison with the classical 
quasilinear equations, the quantum model exhibits a finite-difference structure. The basic 
properties of the quantum quasilinear theory are discussed in Sect. 3, where we derive some 
conservation laws, as well as an appropriate H-theorem for the quasilinear equations. The 
existence of an entropy-like quantity is used to prove that the averaged Wigner function 
relaxes to a plateau, just like in the classical case. However, the distinctive feature of the 
quantum quasilinear equations is the existence of a transient periodic structure in velocity 
space, as shown in Sect. 4. In Sect. 5, the theory is checked against numerical simulations 
of the bump-on-tail and two-stream instabilities, with good agreement with the predictions. 
Conclusions are drawn in Sect. 6. 



II. QUANTUM QUASILINEAR EQUATIONS 

The quasilinear equations for the Wigner-Poisson system were derived long ago by Ve- 
denov 30|, without fully exploring their consequences. For completeness, the derivation 
procedure is reproduced here. The Wigner equation reads 

df df f 

~di + V ~dx = I dV K ^ V ~ V ' x ^^ x ' v '*) ' W 

where 

K(v - v', x, t) = — ^ / d\ e im ^- v '^ n [<j>(x + A/2, t) - <j>(x - A/2, t)] . (2) 
2-kTi J 



Here, f(x,v,t) is the Wigner pseudo-distribution in one spatial dimension, with position x, 
velocity v, and time t. Also, h = h/27r is the scaled Planck's constant, e is the absolute 
value of the electron charge, and m is the electron mass. The electrostatic potential t) 
satisfies the Poisson equation, 

— ( I dvf - n ^) , (3) 



dx 2 e 

where Eq is the vacuum permittivity and n a fixed, neutralizing ionic background. Periodic 
boundary conditions are assumed, with periodicity length L. Accordingly, for any quantity 
A = A(x,v,t), the spatial average is 
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L/2 



(A(x,v,t)) = — dxA(x,v,t). (4) 

L J -L/2 

In particular, it is useful to define F(v,t) = (f(x,v,t)) and to restrict to £)) = 0. 
The quasilinear theory proposes a perturbation solution of the form 

f = F(v,t) + f 1 (x,v,t), <(> = fafat) , (5) 

for small fi and <p±. After averaging the Wigner equation, taking into account that 
(4>x(x,t)) = 0, we obtain from §Q 



OF iem 



dt 2nh 2 



d\dv'e im{v - v ' )x/h 



> L \x- -,t) + 



fi(x,v',t) ) . (6) 



2 ) T V 2 

In the quasilinear theory, the right-hand side of Eq. §6§ is evaluated using the results of the 
linear theory. This implies, in particular, that the mode coupling effects are not taken into 
account. Notice that F changes slowly, since dF/dt is a second order quantity. Thus, for too 
strong damping or instability, the quasilinear theory is no longer valid. Also, the trapping 
effect can be included only in the framework of a fully nonlinear theory. Generally speaking, 
the conditions of validity of the quasilinear theory are still the subject of hot debates at the 
classical level [311 ]. 



We now introduce the spatial Fourier transforms 

1 



f lk (v,t) = -= [ L/2 dxe- ikx h{x,v,t), (7) 

V27T J -L/2 

} lk (t) = -4= [ L/2 dxe'^M^t), (8) 

V27T J -L/2 
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with the corresponding inverse transforms 



fi{x,v,t) 

4>i(x,t) 



k 



(9) 
(10) 



where k = 2im/ L, n = 0, ±1, ±2, ... Linearizing the Wigner equation, Fourier transforming 
it in space and Laplace transforming it in time, we obtain 

/ lfc ( Vj t ) = e ^ lk{t) [ d\ dv' e im ^~ v 'WR( e ifcA/2 _ e -i*A/2) £ ) . (ii) 
2itTi (ujk — kv) J 

In Eq. (tll|), since we are interested only in the long-lived collective oscillations, the initial 
perturbation fi(x, v, 0) was neglected. Finally, u>k stands for the allowable frequency modes, 



obtained from the the well-known 



321 ] quantum dispersion relation 



D(k,iu) = 1 - -£ / dv 



F(v) 







(12) 



n o Jc ~ i^k - kv) 2 - h 2 k A /Am? 
which is assumed to be adiabatically valid, where £ denotes the Landau contour, and uo p = 
(n e 2 /me ) 1 ^ 2 is the electron plasma frequency; for brevity, we omit to write out the second 
argument t of F from here and onward. 

Since f\ and <p\ are real, we have the parity properties 



fi,-k(v, t) = f* k (v, t) , <j> x k (t) = 4>* lk {t) , U_ k = -u 



k ■ 



(13) 



which will be used in the remainder of this paper. By defining 



where Qj~ and 7& are real, the following properties also hold 



Q-k = , 7-Jfc = Ik ■ 



(14) 



(15) 



By using Eqs. (jlip and (TT3]) into ([6]), and carrying out straightforward calculations, we 
obtain 



dF 2nie 2 



dt L 2 h 2 



k^O 



F(v + hk/m) - F(v) F(v -fik/m) - F{v) 
ujk — kv — hk 2 /2m Uk — kv + hk 2 /2m 



(16) 
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Equation (|T6|) can be put in a convenient form by using the parity properties (|T3l) - (|T5l) and 
that |7fc| is small so that jk/ii^k ~ kv ± hk 2 /2m) 2 + 7^] ~ 7vS(Qk — kv ± hk 2 /2m), where 5 
is Dirac's delta function. Hence we express Eq. (|16p as 



dF 
~dt 



An 2 muj 2 ^ i k (t) r 



+ 



n Lh 2 k 2 I 

11 k>0 

hk\ 



he 

2m 



(17) 



in which the growth rate does not appear explicitly, and where we have defined the spectral 
density of the electrostatic field fluctuations as 



e k {t) 



jk 2 \h k (t)\ 2 . 



(18) 



The presence of the delta functions in Eq. (fTTj) emphasizes the fact that only particles 
satisfying the resonance condition Q k — kv ± hk 2 /2m = are taken into account, whereas 
F(v, t) remains unchanged in the non-resonant region of velocity space. 

The time variation of the spectral density is given in the same way as in the Vlasov- 
Poisson case, i.e. 

^ = 2 lk e k . (19) 

Equations (fTTj) and f|T9|) are the quasilinear equations for the Wigner-Poisson system. In 
the formal classical limit % — > 0, they reduce to the well-known quasilinear equations for the 
Vlasov-Poisson system. 



III. PROPERTIES OF THE QUANTUM QUASILINEAR EQUATIONS 

Taking velocity moments of Eq. f[T7|) . several conservation laws can be easily derived. 
For instance, we obtain the conservation of the number of particles, 

^ [ dvF = 0, (20) 



dt „ 

and of the linear momentum 

— / dvmvF = 0, (21) 

dtj K J 

where the dispersion relation f|T2|) and the parity properties of the spectral density have been 
used. In addition, the total energy is also invariant 

d I f , mv 2 2tt v-^ \ , s 

-U dv — F+T ^ it \ =Q . (22) 



Intermediate steps to derive Eq. (122]) require the use of the parity properties as well as the 
second quasilinear equation (1191) . In addition, a useful approximation is to consider fi* 



for k > 0, jointly with 



7fc 



p 



2nnfik 3 



Up hk 
k 2m 



-F 



uj p hk \ 
k 2m J 



(23) 



In order to gain insight into the asymptotic behavior of F(v,t), it is interesting to look 
for an entropy-like quantity. Following Ref. {^J, we consider the quantity J F 2 dv, which, 
from Eq. (fTTj) . can be proved to obey the equation 



ik(t) 

n Lh 2 ^ k 3 



, flu hk\ 



/flu hk 
F | 

k 2m 



< 



(24) 



where the last inequality follows since all terms in the right-hand side are non-positive. This 
result constitutes a sort of H-theorem for the averaged distribution function. 

The time derivative of the non-negative quantity in the left-hand side of Eq. ( 1241) is 
always non-positive, so that asymptotically we have 4 J dvF 2 — > 0, and 



fk \ 
k 2m J 



— F 



flk hk \ 
k 2m J 



0, 



(25) 



for all wavenumbers where the spectral density is non-zero. Equation (1251) also resembles 
the basic equation obtained when applying the Nyquist method to the stability analysis of 
the Wigner-Poisson system 3^]. 

Equation (T2"5"j) is a finite-difference-like version of the classical plateau condition [F'(v) = 
0] in the region where the spectral density is not zero. The finite-difference structure of 
Eq. (125]) favors the appearance of oscillations in velocity space, which are not present in the 
classical case. These oscillations are analyzed in detail in the next section. 



IV. TRANSIENT QUANTUM OSCILLATIONS IN VELOCITY SPACE 

In the derivation at the end of Sect. 3, it was implicitly assumed that the spectral density 
is not zero in a broad region of momentum space. In contrast, let us see what happens in 
the idealized situation where the spectral density is strongly peaked at a single mode 
K. Denoting the associated spectral density by e(t) and using the growth rate (1231) . the 
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quantum quasilinear equations read 



dF mLu; 



hK 2 



and 



di(t) Trmu 



u p hK\ ( u p HK\ 



K 2m) \K 2m J 



(27) 



dt riohK 3 

where, for simplicity, the approximation Qk = u p was also adopted. 

A particular class of stationary solutions of Eqs. (I26i)-(l27l) is given by any function F(v) 
that is a periodic in velocity space, with period hK/m: 

F(v + hK/2m) — F(v — TiK/2m) = 0. (28) 

Assuming F(v) ~ exp(iav) in Eq. (1281) . with a to be determined, we easily obtain the 
characteristic equation sin(a^i^/2m) = 0. Hence, the general (exact) equilibrium solution 
is the linear combination 

f 2TTnv\ / 2Tcnv\ 

F(v) = a + 2_^a n cos I — - — 1 + > y 6 n sin I — - — 1 , (29) 

?1=1 ^ u / n= l \ « / 

where arbitrary real constants. Notice the singular character of the quantum oscil- 

lations, whose "wavelength" of the fundamental mode (n = 1) in velocity space, X v = hK/m, 
tends to zero as h — > 0. The solution given by Eq. (1291) represents periodic oscillations in 
velocity space. However, this is necessarily a transient solution that cannot be sustained 
for long times. Indeed, strictly speaking, Eq. f[2"9l applies only at the resonance, which is 
a set of measure zero on the velocity axis. The generation of harmonics with wavenumbers 
k 7^ K (which is forbidden in the quasilinear theory) would ultimately lead to a broad energy 
spectrum. 

When many wavenumbers are present, the simple periodic solution given in Eq. (1281) does 
not hold anymore, because different values of k induce different velocity-space wavelengths 
A„. Only the solution F(v) = const, holds independently of k. Therefore, we expect that a 
plateau will eventually appear on a finite region in velocity space. Nevertheless, Eq. ( 1291) 
provides an estimate for the characteristic oscillation length of the averaged Wigner function 
in velocity space, near resonance. 
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FIG. 1: Initial velocity distribution used in the simulations of the bump-on-tail instability, both 
for the Vlasov and Wigner cases. 

The above arguments also suggest that monochromatic waves are the best candidates for 
displaying such quantum oscillations. In contrast, if the energy spectrum is broad from the 
very start, the formation of a plateau is likely to be favored over the appearance of peri- 
odic oscillations. In the following section, these predictions will be compared to numerical 
simulations of the Wigner- Poisson system. 

V. NUMERICAL SIMULATIONS 

We have performed numerical simulations of the Vlasov and Wigner equations using a 
phase-space code based on a splitting method |35j . In the Wigner equation, the acceleration 
term (jSJ) is a convolution product in velocity space and is therefore calculated numerically 
by Fourier transforming it in velocity space. 

To study the differences in the nonlinear evolution of the Wigner and Vlasov equations, 
we have simulated the well-known bump-on-tail instability, whereby a high-velocity beam 
is used to destabilize a Maxwellian equilibrium. We use the initial condition / = (1 + 
<5)(no/v / 27Tfth)[0.8 exp(— v 2 /2vj h ) + 0.4exp(— 2(v — 2.5vth) 2 /v^ h )}, where 5 represents random 
fluctuations of order 10~ 5 that help seed the instability (see Fig. 1). Here v t h = \/k B T e /m 
is the electron thermal speed. We use periodic boundary conditions with spatial period 
L = AOnXoe, where Xn e = v t h/uj p is the Debye length. Three simulations were performed, 
with different values of the normalized Planck constant, defined as H = Tioj p /mvl h : H = 
(Vlasov), H = 1, and H = 2. 
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H = C H=l H=2 

5.0 1 5.0 1 5.0 




FIG. 2: Simulations of the Wigner-Poisson and Vlasov-Poisson systems, at time co p t = 200, for H = 
(Vlasov, left frame), H = 1 (middle frame), and H = 2 (right frame). Initially monochromatic 
spectrum. Top panels: electron distribution function f(x,v) in phase space. Bottom panels: 
spatially averaged electron distribution function F(v) in velocity space. 

In order to highlight the transient oscillations in velocity space, we first perturb the 
above equilibrium with a monochromatic wave having Xn e k = 0.25 (i.e., a wavelength of 
87rA.D e ). Figures 2 shows the results from simulations of the Wigner-Poisson and Vlasov- 
Poisson systems. In both simulations, due to the bump-on-tail instability, electrostatic 
waves develop nonlinearly and create periodic trapped-particle islands (electron holes) with 
the wavenumber k = 0.25. The theory described in the previous sections predicts the 
formation of velocity-space oscillations in the Wigner evolution, which should be absent in 
the classical (Vlasov) simulations. This is the case in the results presented in Fig. 2, where 
the oscillations are clearly visible. 

In order to estimate their wavelength, a zoom on the spatially averaged electron distri- 
bution function F(v) is shown in Fig. 3. According to the quasilinear theory, the velocity 
wavelength should be equal to A„ = hk/m for the fundamental mode with n — 1 [see Eq. 
(l29l)]. In our units, this yields X v /vth = HkXoe, which is equal to 0.25 for H = 1 and to 0.5 
for H = 2. The wavelengths observed in the simulations are slightly smaller: X v /v t h — 0.17 
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v / v tl 




v / v tl 



FIG. 3: Zoom on the spatially averaged electron distribution function F(v) in velocity space for 
the same cases as in Fig. 2. 

and 0.35 for H = 1 and H = 2, respectively. However: (i) the oscillations are absent in 
the Vlasov expected, (ii) the order of magnitude of the wavelength is correct, and 

(hi) the wavelength is proportional to H, in accordance with the quasilinear theory. The 
slight discrepancy in the observed value of A„ may have at least two origins. First, spatial 
wavenumbers different from 0.25 can be excited due to the nonlinear mode coupling. Second, 
other modes with n > 1 [see Eq. (1291) ] can affect the velocity-space wavelength. 

When the initial excitation is broad-band (i.e., wavenumbers 0.05 < k < 0.5 are excited), 
ihe electron holes start merging together at later times due to the sideband instability 



36 



371 ] (see Fig. 4). At this stage, mode coupling becomes important and quasilinear 



theory is not capable of describing these effects. As the system evolves toward larger spatial 
wavelength, the evolution becomes progressively more classical, with the appearance of a 
plateau in the resonant region. Nevertheless, at u p t = 500 the Wigner solution still displays 
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t = 500. t = 500. 




FIG. 4: Simulations of the Wigner-Poisson and Vlasov-Poisson systems, for u p t = 500, for H = 
(Vlasov, left panel) and H = 1 (Wigner, right panel). Initially broad wavenumber spectrum. Top 
panels: electron distribution function f(x,v) in phase space. Bottom panels: spatially averaged 
electron distribution function F(v) in velocity space. 

some oscillatory behavior in velocity space, which is absent in the Vlasov evolution. 

Another set of simulations were performed for the case of a two-stream instability. The 
initial distribution function is composed of two Maxwellians with thermal speed v t h, each 
centered at v — ±2vth (see Fig. 5). Only the fundamental mode of the system, with the 
wavenumber K = 0.2\]^ e , is excited, and it grows exponentially due to instability. Three 
simulations were performed, with different values of the normalized Planck constant: H = 
(Vlasov), H = 1, and H = 2. The expected oscillation wavelength in velocity space is 
A„ = TiK/m. For the three simulated cases, it should take the values A^ = 0, A„ = 0.2v t h, 
and A^ = 0Av t h- A zoom of the averaged Wigner function F(v, t) around the velocity v — 
is plotted in Fig. 6 for the three cases, at time uj p t = 110. The velocity-space oscillations 
are absent from the Vlasov simulation. In the Wigner cases, their wavelength is rather close 
to the theoretical value; in particular, it appears to grow with the scaled Planck constant, 
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FIG. 5: Initial velocity distribution used in the simulations of the two-stream instability, both for 
the Vlasov and Wigner cases. 

as expected from the theory. As in the bump-on-tail case, the oscillations tend to disappear 
over longer times. 



VI. CONCLUSION 



In this paper, the quasilinear theory for the Wigner-Poisson system was revisited. Con- 
servation laws and the asymptotic solution were established. Distinctive quantum effects 
appear in the form of a transient oscillatory behavior in velocity space. Such quantum ef- 
fects are favored when the energy spectrum is restricted to few modes. For longer times, the 
plasma tends to become classical - due to the spatial harmonic generation and the mode 
couplings - at least as far as the averaged Wigner function F(v,t) is concerned. Thus, 
just as in the classical case, F(v , t) evolves asymptotically toward a plateau in the resonant 
region of velocity space. It would be interesting to investigate whether monochromaticity 
enhances quantum effects in plasmas in general, which would represent a useful property for 
an experimental validation of the quantum plasma models. 
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